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The statistical properties of interstellar turbulence are studied by means of three- 
dimensional high-resolution HD and MHD simulations of a SN-driven ISM. It is 
found that the longitudinal and transverse turbulent length scales have time aver- 
aged (over a period of 50 Myr) ratios of 0.5-0.6, almost similar to the one expected 
for isotropic homogeneous turbulence. The mean characteristic size of the larger 
eddies is found to be ~ 75 pc. Furthermore, the scalings of the structure functions 
measured in the simulated disk show unambiguous departure from the Kolmogorv 
(1941) model being consistent with the latest intermittency studies of supersonic 
turbulence (Politano & Pouquet 1995; Boldyrev 2002). Our results are indepen- 
dent of the resolution, indicating that convergence has been reached, and that the 
unresolved smaller dissipative scales do not feed back on the larger ones. 

1. Introduction 

Interstellar turbulence is mainly driven by the energy injected into the ISM 
by supernovae, with the driving scale still being uncertain. It is also unclear 
what the statistical properties of the turbulent interstellar gas are, if the full 
available range of energies is taken into account. In this paper we discuss 
the statistical properties of the interstellar turbulence in the Galactic disk 
based on three-dimensional adaptive mesh refinement (AMR) simulations 
of the ISM, which include the disk-halo-disk circulation In particular 
we explore the injection scales (section 2) and the scalings of the velocity 
structure functions (section 3) of the interstellar turbulent gas and discuss 
(section 4) their implications. 

2. The injection scale of interstellar turbulence 

The outer scale of the turbulent flow in the ISM is related to the scale at 
which the energy in blast waves is transferred to the interstellar gas. Such 
a scale can be determined by means of the longitudinal and transversal 
correlation lengths, L\\ and Lkk, respectively. Here, k = 2,3 refer to 
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the directions perpendicular to the 1-direction along which the correlation 
lengths are calculated. For isotropic turbulence Lkk — 0.5Ln. Figure 1 
shows the history of in (left panel) and ^22/^11 (right panel) during 50 
Myr of evolution of the unmagnctized and magnetized ISM. For details 
on how Ln and L 22 are calculated see [ 6 ]. Although (I u ) ( <~ 75 pc, the 
scatter of Ln around its mean results from the oscillations in the local 
star formation rate, where the formation and merging of superbubbles, is 
responsible for the peaks observed in the two plots. The similarity between 
the average HD and MHD injection scales is due to the fact that magnetic 
pressure and tension forces cannot prevent break-up as long as L < (3pX, 
where L and A are the scale lengths of thermal and magnetic pressures 
(including tension forces) and (3p = AnP/B 2 is the plasma beta. 



- MHD: <L M > = 75.2pc 

- HD: <L n > = 72.7pc 
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_ MHD: <L 22 /L n > = 0.6 
- HD: <L 22 /L n > = 0.Sl ~ 




Figure 1. History of the characteristic size (given by Ln) of the larger eddies (left 
panel) and of the ratio Z/22/L11 (right panel) for the HD (dashed line) and MHD (solid 
line) runs. 



Despite the large scatter seen in L22/L11, the time average {L22/ L\\)t 
over the 50 Myr period is 0.51 and 0.6 for the HD and MHD runs, re- 
spectively. The discrepancy from 0.5 by about 20% in the MHD case is a 
consequence of the anisotropy introduced by the field into the flow. The 
(L 22 /Lu)t ~ 0.5 in the HD run indicates that in a statistical sense the 
interstellar unmagnetized turbulence is roughly isotropic. 

3. Scalings of the Structure Functions 

The statistics of turbulent flows in physical space is commonly characterized 
by the velocity structure functions of order p defined as S p (l) = (\AVi\ p ), 
with AVJ = v(x + l) —v(x), where v(x + l) and v(x) are the velocities along 
the x— axis at two points separated by a distance I, such rj <C ! -C L, L be- 
ing the energy integral scale and r\ the Kolmogorov scale, respectively. For 
homogeneous and isotropic turbulence, the Kolmogorov (K41) [ 2 ] theory 
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predicts that (| AX^| 3 ) = —\d for I ^ t/ 3 / 4 e -1 / 4 . Here, v is the kinematic 
viscosity and ( ) stands for average over the probability density function of 
AV(l). As a consequence of this relation, within the inertial range one can 
write (|A^n oc (|A^| 3 ) c(p) . Experiments show that this expression, with 
the same scaling exponents, is valid in a wide range of length scales for large 
as well as small Reynolds numbers, even if no inertial range is established 
[ 7 ]. Figure 2 displays, for the unmagnetized interstellar gas, (| AV| } (top 
of left panel) and {| AT/| 6 ) (bottom of left panel) as functions of (| Al^l 3 ) 
and its best fits in dashed lines with slopes 0.73 and 1.57 corresponding 
to the scalings of £(2) and C(6), respectively. These scalings arc different 
than those predicted by the K41 theory given by p/3 and similar to those 
proposed by Boldyrev [ 4 ]. More generally, the structure functions can be 




Figure 2. Left: Log-log plot of (|AV| 2 ) (top panel) and (|AV| 6 ) (bottom panel) as 
function of (|AV| 3 ) and its best fits in dashed lines. Right: Comparison between theory 
(plus MHD simulations) and experiments of the exponents £(p) vs. order p. 



written as S p (l) oc l^ p \ where ((p) is given by Cb)/C(3) = 1P+C (1 - TP), 
with C = (1 — 37)/ (l — S 3 ) and C(3) = 1; £ is a measure of the degree of 
intermittency of the flow, C (the Hausdorff codimension of the support of 
the most singular dissipative structures) is related to the structures' Haus- 
dorff dimension D by C = 3-D. In incompressible flows this corresponds to 
one-dimensional vortex filaments [ 3 ], hence C — 2, while the most singular 
dissipative structures in supersonic turbulence are two-dimensional shock 
fronts [ 4 ], and in fully developed MHD turbulence these are 2D micro- 
current sheets [ ]; hence C = 1. 7 is related to very high order moments 
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and has the value of 1/9 for incompressible flows [ 3 ], with a similar value 
quoted by [ 4,s ] for supersonic HD and MHD turbulence. Hence, S 3 = 2/3 
or 1/3 if C = 1 or 2. The degree of intermittency in the compressible 
turbulence is higher than in the incompressible case. In Kolmogorov tur- 
bulence £(p) — p/3, which implies that C = independently of the value 
of S. A comparison between the different theoretical scalings and the ones 
derived from the MHD run is shown in Figure 2. Although the relative 
scalings, which express the dependence of the moments at different orders, 
are universal, they show unambiguous departure from the Kolmogorov law 
due to intermittency of interstellar turbulence. For further discussion on 
these results see [ 6 ]. 

4. Final Remarks 

In subsonic turbulence energy is injected at the outer scales and transferred 
to the smallest scales, with dissipation by molecular viscosity setting in at 
the Kolmogorov inner scale. Dissipation is a passive process as it proceeds 
at a rate determined by the inviscid inertial behaviour of the large eddies. 
Such an energy cascading corresponds to a divergence-free behaviour of 
the flow in the inertial range. However, the ISM is a compressible medium 
swept up by shocks, which are in general more efficient in dissipating energy, 
although most of it occurs again at small scales through low Mach number 
shocks. 
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